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Abstract 



^ ' We present a calculation of the cross section and the event generator of the 

^^ ■ reaction vre — > 7re7. This reaction is sensitive to the pion generalized dipole 

Q ■ polarizabilities, namely, the longitudinal electric aL{q^), the transverse elec- 

C^^ , trie axiq'^), and the magnetic I3{q'^) which, in the real-photon limit, reduce to 

,__! I the ordinary electric and magnetic polarizabilities a and f3, respectively. The 

^D ' calculation of the cross section is done in the framework of chiral perturbation 

,J^ ■ theory at 0{p'^). A pion VCS event generator has been written which is ready 
for implementation in GEANT simulation codes or for independent use. 
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I. INTRODUCTION 

Pion electromagnetic polarizabilities are important dynamical properties, providing a 
window to the pion's internal structure (for an overview see, e.g., Refs. |jl|). In a classical 
picture, polarizabilities determine the global, i.e. integral, deformation response of a com- 
posite system in static, uniform external electric and magnetic fields. As first shown in Refs. 

0, a determination of hadron polarizabilities is possible via 7-hadron Compton scattering. 
So far, pion electric (a) and magnetic (/?) polarizabilities were investigated via 77r -^ 771 
real Compton scattering (RCS). Since pion targets are unavailable, pion RCS was done 
using different artifices, such as high-energy pion- nucleus bremsstrahlung, 71^ Z ^ ti^Z'j 

01, radiative pion photoproduction off the proton 7p -^ ^n^n Q, and the crossed-channel 
two-photon reaction 77 —>■ vrvr as embedded in the e~e^ —>■ e^e^TT~^TT~ process j^,^] (see Fig. 
|l]). The values of a extracted from these experiments are given in Table |. 

As can be seen, the results scatter substantially and have large uncertainties, and there- 
fore, more experimental efforts are evidently needed. A new radiative pion photoproduction 
experiment at the Mainz Microtron MAMI has recently taken data [0 , another one has been 
proposed at Jefferson Lab |p, an experiment investigating the inverse reaction ir^p -^ wyy 
has been proposed at TRIUMF 0], a new 77 -^ ini experiment was proposed at Frascati 
PI| , and new Primakoff (bremsstrahlung) experiments were proposed at CERN COMPASS 
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Electromagnetic studies with virtual photons are advantageous, since the energy and 
the modulus of the three-momentum of the virtual photon can be varied independently. 
Furthermore, one also gets access to longitudinal degrees of freedom which cannot be tested 
by real photons. In the pion VCS reaction 7*7? ^ 771 with a virtual space-like initial-state 
photon and a real final-state photon, the so-called generalized polarizabilities (GPs) can be 
measured in the space-like region. Originally, such quantities were introduced in the context 



of nuclei |]T2| and have recently been reinvestigated for the nucleon in Ref. [I^. The analysis 
of the nucleon GPs is based on a multipole expansion in the center-of-mass frame below 
the pion production threshold. The first experimental results for generalized polarizabilities 
of the proton have been reported by the VCS collaboration at MAMI WM and are in good 
agreement with the predictions of chiral perturbation theory ||15[ . 

Recently, another approach to the analysis of the VCS tensor has been performed in a 
completely covariant framework |ll6| , p!7| . In particular, the reaction 7*7r — > jn is described in 



terms of three Lorentz-invariant functions, each depending on three scalar variables. In the 
limit of a soft final real photon, q' — > 0, the structure-dependent response is encoded in three 
generalized dipole polarizabilities ^^(g^), ^^(g^), and /3(g^) [|l7l. The Fourier transforms of 
these functions are related to the electric polarization and magnetization induced by soft 
external fields fl^. In other words, the VCS generalized dipole polarizabilities allow one to 
access the local polarization response, whereas the RCS polarizabilities a = 01,(0) = ariO) 
and P = /3(0) only provide the integrated, i.e. global, response. Furthermore, two electric 
polarizabilities are needed to fully reconstruct the electric polarization, the longitudinal 
electric polarizability being related to the divergence of the induced electric polarization, and 
the transverse electric polarizability describing rotational displacements of charges which do 
not result in a modification of the charge density (such as, e.g., a rotating spherical charge 
distribution). On the other hand, due to the transverse character of the magnetic induction. 



only one generalized magnetic dipole polarizability appears. 



Motivated by the success of the nucleon VCS experiment ||1J] , by studying the reaction 



7re 



7re7, 



(1) 



we attempt to get a hold of the pion VCS matrix element 7*(g) + vr(pj) — > 7(g') + 7r(p/), 
which can then be utilized as a new experiment for a pion polarizability measurement. 

For the numerical analysis of the above reaction, we make use of the results obtained in 
the framework of chiral perturbation theory at 0{p^) [l^. At this order in the momentum 
expansion, the generalized dipole polarizabilities display a degeneracy, /?(g^) = —aLil'^) = 
— ar('?^) which will be lifted at 0{p^) and 0{p^), respectively. 

Our work is organized as follows. In Sec. II, we briefly discuss the predictions for the 
generalized dipole polarizabilities as obtained in the framework of ChPT at 0{p'^) [0. We 



provide a physical interpretation of the longitudinal generalized electric dipole polarizability 
in terms of the induced polarization charge [|17i. Section III briefly describes the kinematics 



and notation used. Section IV contains a discussion of the various pieces entering the cross 
section, whereas some general features of the cross section are discussed in Sec. V. Section 
VI introduces the ingredients of the event generator. In Sec. VII we discuss the sensitivity of 
the cross section to the polarizabilities. Section VIII contains the results for the cross section 
obtained after Monte Carlo integration. The present experimental status is summarized in 
Sec. IX. In Sec. X we summarize our results and draw some conclusions. The analytical 
results for the squared invariant amplitudes are given in the Appendix. 

II. VIRTUAL COMPTON SCATTERING AND GENERALIZED DIPOLE 

POLARIZABILITIES 



A new frontier for polarizability studies is the subject of VCS |]T^,|T8|,|T9[ . There are two 
different ways of studying this reaction in different kinematical regimes regarding the virtual 
photon. In the first case, the initial photon is real and the final virtual photon fragments into 
This corresponds to positive g'^ > 4ml ^^"^ "^ill '^'^^ be discussed here, 
we concentrate on the second case, where the initial virtual photon is produced 



a Dalitz pair pO 
Rather 



in a particle scattering process, thus leading to q < 0, and then scatters from a target 
to a real final-state photon. This process probes pion structure via the generalized dipole 
polarizabilities aiiq"^), ariq"^), and Piiq'^) which in the real-photon limit, g^ — *> 0, reduce to 
the RCS polarizabilities ai^O) = «t(0) = a and /5(0) = /? [p!6| , p^ . The predictions in the 
framework of ChPT at O(p^) read |T^ (see Fig. |^) 
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The predictions for the charged pion RCS polarizabihties result from an old current- algebra 



low-energy theorem |21 
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where the low-energy constants are fixed via the radiative pion decay vr"*" —>■ e^Vel- Correc- 
tions at 0{p^) were shown to be 12% and 24% of the 0{p^) values for a and /5, respectively 
22 1 . The q^ evolution of Eq. (^ [and also of (|^)] is completely determined in terms of the 
pion mass m.,^ and the pion-decay constant F^^ = 92.4 MeV. 

The generalized dipole polarizability aL(5'^) can be interpreted in terms of the influence 



of a static, uniform electric field E on the pion charge form factor |jT^. In the rest of this 
section we will not make use of a covariant notation. Breit frame kinematics has to be applied 
in the transition to the relativistic case. The presence of a static uniform electric field E will 
induce an electric polarization, the divergence of which is defined as the polarization charge 

density. 
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The potential created by the original charge distribution in combination with the induced 
electric polarization will look like the one which is due to an effective charge density: 



Peff(^) = po(^) - V ■ P{x) 



(7) 



where po{x) is the charge density without external electric field. The form factor of this 
effective charge density is then given by 



F{q) = Jd^xexp {-iq- x)p,s{x) = Fo{q) + AF{q), 



(8) 



where -Fo(^ is the form factor of the "undistorted" charge distribution, and 

AF{q) = Amq- EaL{q) (9) 

is the lowest-order change in the form factor due to the external electric field. 

In the above discussion we have concentrated on the polarization charge density, i.e., the 
divergence of the electric polarization, or, in terms of its Fourier transform, the longitudinal 
contribution, only. A full recovery of the polarization requires the transverse piece as well. 
For details, the interested reader is referred to the comprehensive discussion of Ref . [^ . 

III. VCS KINEMATICS AND DIFFERENTIAL CROSS SECTION 



One possibility of studying the pion VCS matrix element would consist of investigating 
the electron scattering reaction ep -^ e''jTT^n as an extension of the RCS experiment 7p ^ 



7vr"'"n 0J3]- ^^ s^*^^ ^ reaction, the VCS contribution is hidden as a subprocess in terms of 
a pion-pole diagram in the t channel. However, the initial pion is off its mass shell and thus 
the extraction involves an extrapolation to the unphysical region with all its ambiguities. 

In the present work, we consider the pion VCS reaction Ti^Pi) + e{k) — > n{pf) + e{k') + 
7(g', e') which, in principle, has the advantage of a direct access to the VCS matrix element 
without any extrapolation to the pion pole. Here the pion four-momenta in the initial and 
final states are denoted by pf = (_Ej,pj) and pj = {Ef,'pf), respectively. Similarly, the 

electron four-momentum reads k^ = {Ek,k) [k'^ = {Ek',k')] and the final-state photon is 
characterized by its four-momentum g'^ = {Eqi,q') and polarization vector e''^ = (e'°,e*'). 
The starting point of our discussion is the differential cross section in the convention of Ref. 

d(7 



m: 



jEk'Eq, ^(p. . fc)2 _ ^2^2 {27ry 



X \M\^\p^ + k-pf-k' - q')d^pfSk'd^q', (10) 

where the invariant amplitude M. contains the complete information on the dynamics of the 
process. As usual, the quantity |A1p indicates an average over the initial polarizations and 
a summation over the final ones. 

The cross section is more conveniently expressed in terms of the five independent invariant 
variables (see Fig. ^ 

s= {pi + kf, 
si = {k' + q')\ 
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This introduces a Gram determinant [24 
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into the expression for the fourfold differential cross section 
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The physical region is determined by the condition A4 < 0. In an actual experimental 
realization, the energy of the pion beam is fixed such that s = 7nl + 7nl + 2Ei7n^ in terms 
of laboratory variables. In other words, the differential cross section of Eq. (|13D then only 
depends on four variables. 



IV. CALCULATION OF THE SQUARED INVARIANT AMPLITUDE 

At lowest order in the electromagnetic coupling, the total invariant amplitude of the 
reaction 7r^(pj) + e~{k) -^ 'K^{pf) + e~{k') + 7(g') can be divided into Bethe-Heitler (BH) 
and VCS contributions, M. = A^bh + -^vcs (see Fig. ^). In the case of the BH diagrams of 
(a) and (b), the real photon is emitted by the final and initial electrons, respectively. The 
hadronic part of the corresponding amplitude is completely determined in terms of the pion 
electromagnetic form factor. We denote the space-like four-momentum of the virtual photon 
emitted by the pion transition current in the BH process hj r = pi—pf = q' + k' — k. Due to 
the hermiticity of the electromagnetic current operator, the electromagnetic form factor is 
real in the space-like region [^]. Applying the Feynman rules of QED, it is straightforward 
to evaluate the BH amplitude. 

The pion VCS amplitude is part of the VCS diagram (c) of Fig. |^. It is this contribution 
which contains the structure information in terms of the generalized dipole polarizabilities 
(and, of course, also of higher-order contributions). 



According to Low's theorem p6|, the irregular terms of the amplitudes A^bh and A^vcs 
behave as 1/Eqi, where the singularities originate from the emission of the (soft) final photon 
at external lines. Gauge invariance then fixes terms of (9(-E'°,) which do not depend on the 
direction q' (for a detailed discussion, see Refs. P7|P^ ). Model-dependent terms in Alycs 
are at least of order Eg'. In this context, the VCS amphtude is divided into generalized 
Born terms A^Born, containing ground-state properties in terms of the electromagnetic form 
factor, and a residual non-Born contribution A^nb, describing model-dependent internal 



structure (see Refs. |ll6|Jl7| , |28|| ) . In particular, the separation is performed such that both 



pieces are separately gauge invariant and the generalized Born terms contain all soft-photon 
singularities. The non-Born part is of 0{Egr). 

When squaring the total amplitude, the irregular terms of |Aip are of order E~,'^ and 
E~^. Such terms entirely originate from the BH part and the generalized Born terms of the 
VCS part. In particular, these contributions are completely known and can be expressed 
in terms of the electromagnetic form factor. Recall that the space-like virtual-photon four- 
momenta are r and q in the Bethe-Heitler and VCS contributions, respectively, and only in 
the limit q' -^ they coincide. 

The spin-averaged squared invariant amplitude can be decomposed into three parts, a 
pure VCS part, a pure BH part, and the interference term between BH and VCS: 

1 



\M\' = 7; E(-^vcs + A1bh)(A<^cs + -^bh) 



= IA^bhP + lAlvcsP + A^vcsAIbh + A^bhA^vcs- (15) 

The Lorentz invariant matrix element for the BH amplitude reads 



le^ 



Mbu = ±^e'*^{q', A') [{Pf + p.YF (r^)] L^,^ (16) 

where the upper (lower) sign applies in the case of a 7r+ (tt^) beam. Here, L^^ is nothing 
else but the second-rank electron Compton tensor. 
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Equation (|T^ nicely displays the origin of the singularity for g' — > as being due to the 
internal electron lines of the electron Compton scattering s- and w-channel diagrams ap- 
proaching the on-shell limit. The squared matrix element decomposes into the product of a 
pure leptonic and a pure hadronic part, having the same overall sign for both vr"*" and vr", 
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Since l^^ is a product of a term containing the final-state photon polarization vector times 
another term containing its complex conjugate, the sum over the polarizations can be carried 
out immediately yielding the well-known result 



E^;*(^')<(A') 
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(20) 



This simple version holds true only for the case of conserved quantities, i.e., quantities 
which vanish after replacement of the polarization vector with the corresponding photon 
momentum.^ This is the case in the present situation. Contracting the electron spinors and 
using trace techniques according to Ref. p3[, l^^ can be expressed in terms of the four- 
momenta of the initial and final electrons k and k', respectively, and the final-state photon 
four- momentum q'. Finally, using on-shell conditions for the electron, k"^ = k' "^ = ml, we 
obtain 



/BH _ _ 

"•^ 2mlk ■ q 



— -^ {2k^k^ - q'^kl - k'^ql + h,,^ + k ■ k'g^,^ 



2m1k' ■ q' 



- - \2k^k^ + q^ky + fc^g^ + /i^jy + k ■ k g^y 



2m1k ■ q'k' ■ q' 



k ■ k' [q'^kl + k'^ql - q'^k^ - k^ql - 2h^y - 2mlg^y) + 2mlq'^ql 



^For the general case, one would have to to use the complete formula given, e.g., in Ref. |2S] 
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(22) 



Note that l^^ is symmetric under k ^->- —k' as well as under fi ^->- u. Also r'^l^^ = = r'^l^^- 
The analytical result for the contraction of Eq. (|18D is given in Eq. (|A2|) of the Appendix. 

The VCS amplitude is obtained by contracting the pion VCS tensor with the initial 
virtual-photon and final real-photon polarization vectors e^ = eu'')^u/q^ and e'* , respectively, 



le 



-Mvcs = —tU {k', s') ^,u {k, s) e';^rcs- 



q 



(23) 



The generalized Born terms of the pion VCS tensor contain, in addition to the s-channel 
and u-channel pole terms, a "contact term" in order to establish gauge invariance |T^,[T^J2^ 






{2p, + qr{2pf + qr ^ {2pf - qY(2p,- q') 



S2 -m^ 



U2 -m^ 
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(24) 



where U2 = {pi — q'Y- The most general form of the non-Born tensor involves three inde- 
pendent structures with corresponding invariant amplitudes depending on three invariant 
variables (see Refs. |T6| , p!7| , |30| , |3T| ) . Here we will only consider the low-energy limit and make 
use of the 0{p'^) predictions of Eq. @). In this approximation the non-Born terms read 



M^'' 



iTim 



NB 



-(gV-g-gy^)«L(g'), 



(25) 



where aii^q^) is given in Eq. (^. 

The squared pure VCS matrix element involves the standard lepton tensor known from 
the one-photon exchange approximation in electro-production processes, 

%^=2 I^t""*^^'' ^'hM^^ s)] H^'^ shuu{k, s)]* 



2ml 

1 
4m? 



k^kjy + /c^/Cy + g^,y[m^ k ■ k ^ 



(26) 



where K = k + k' and q = k — k'. The second equation is particularly useful when contracting 
with another "conserved" tensor. The lepton tensor is contracted with a hadronic tensor of 
the form0 



^At C(p^), for space-like q'^ all quantities are real, even for S2 > (Smj^)'^, and we could, in principle, 
omit the complex conjugation. However, for the sake of generality, we stick to the full notation. 
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^rcs = E ^';MZ:sKM'(rcsr = -^P.-Mrcs^vcs, (27) 
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where, again, we made use of Eq. (^). Using Eq. (|26|) together with current conservation, 
we then obtain 
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Squaring the VCS tensor as in Eq. (pSf ) gives rise to three distinct contributions. They 
involve, respectively, the contraction of the generalized Born terms, of the interference be- 
tween the generalized Born terms and the residual terms, and of the residual terms involving 
the polarizabilities. The respective analytical results can be found in Eqs. (|A^), (^^), and 
( |X5| ) of the Appendix. 

The interference term between the BH and the VCS amplitude can, again, be written 
as the contraction of a hadronic and a leptonic part. Since the BH amplitude changes sign 
under the substitution vr^ ^^ vr", the interference term is sensitive to the pion charge. The 
upper and lower signs apply to the cases of a vr"*" and a tt~ beam, respectively, 

1 
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L^^^, = lj:^{k,shMk',s')Ll^. 



Using trace techniques one obtains 
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1 fC r 

+gvp[ijnl - k- k')q'^ + k' ■ q'k^ - k ■ q'k'^ \ ^Tr~, ^''K + ^p^^ + ("^e - k ■ k')g^ 
-^Jk^-k']. (30) 
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Finally, the interference term can be divided into the interference of the BH part with the 
VCS Born and VCS non-Born contributions, respectively. The results are given in Eqs. (|A7|) 
and (|A8|) of the Appendix. 



V. GENERAL FEATURES OF THE VCS DIFFERENTIAL CROSS SECTION 

Some general features of the fourfold differential cross section, Eq. (|T3p, can be inferred 
from Eqs. ( |A1| ) - (^) of the Appendix. The Si dependence is dominated by the (si —ml)~^ 
pole of the BH amplitude and the cross section varies according to this law with only a 
slight modification through the si dependence of the VCS amplitude. The S2 dependence is 
dominated by the (s2 — rri^)'^ pole of the VCS amplitude with modifications due to the S2 
dependence of the BH amplitude, but in this case the modification is not so small as in the 
case of Si. 

The r^ dependence is determined by the l/(r^)^ pole of the squared BH amplitude, and 
the g^ dependence follows the l/(g^)^ pole behavior typical of electron scattering. The energy 
of the outgoing pion is expected to be high and the angle is expected to be small according 
to l/(r^)^ behavior of the cross section. The energy of the final-state photon is mainly 
expected to be low, as a consequence of the infrared divergence of the BH amplitude. The 
angular behavior of the final photon is determined by the l/(si — ml) and l/(s2 — wi^) poles 
of the BH and VCS amplitudes, respectively. The more interesting photons related to the 
GPs are expected to have higher energies. The behavior of the final electron is completely 
determined by the l/(g^)^ behavior of the cross section. 

The different behavior under the substitution vr^ — > 7r+ of A^bh and A^vcS; 

-MbhIvt") = -^BH(7r+), 

-Mvcs(vr-)=A^vcs(vr+), (31) 

may be of use in identifying the contributions due to internal structure by comparing the 
reactions involving tt~ and tt~^ beams for the same kinematics [|T9|] : 



dcx{n+) - da{n-) ~ 4i?e{A<BH(vr+)A^^cs(^+)}- (32) 

This circumstance suggests the possibility of using such a subtraction to remove from the 
data sample the pure BH and VCS contributions. 



VI. VCS EVENT GENERATOR 
The VCS event generator has been written based on the differential cross section of Eq. 



(O) in combination with the matrix elements of Eqs. (|A2|)-([X8D, using three-body- final-state 



kinematics. The acceptance- rejection method [32| is used for event generation 



A. Generation of the invariants 

The generation starts with the calculation of the squared invariant mass s associated 
with the incoming pion beam and the electron target [see Eq. (]TT|)]. In terms of laboratory 
(lab) variables, s is given by 

s = ml + ml + 2Eime, (33) 
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where Ei is the incident beam energy. 

The VCS cross section increases rapidly if either the direction of the outgoing real photon 
is close to the direction of one of the outgoing charged particles (due to the l/(si — ml) 
and l/(s2 — ml) poles in the BH and VCS amplitudes, respectively), or if the energy of the 
outgoing real photon approaches zero (due to the infrared divergence of the l/(si — ml — 
q^ + r'^) BH pole). 

Therefore, if events are generated in the pole region, the efficiency of the acceptance- 
rejection method for the more interesting non-pole regions can be rather low. In order to 
generate interesting events at an acceptable rate, the pole regions are cut. Invariants are 
generated in the following regions: 

1000 m^ < si< mj, 

2ml < S2< m^, 

ql,, < g2< 2ml - '2m,E^'^, 

rLn < r'< -'^rn.Ef^ + g' + Si - ml (34) 

We choose -E"™'" = 10 GeV in order to cut the 1/g^ pole and to understand the trigger 
acceptance. For the minimal photon energy we choose E'T" = 5 GeV in order to cut the 
infrared divergence of the BH amplitude. For the generated invariants si, S2, q^-, and r^, 
the positivity of — A4 of Eq. ([T^) is checked. In case of a positive A4, the generated event 
is rejected and a new event is generated. Actually, the requirement of a negative A4 is 
equivalent to the requirement of energy- momentum conservation. For an accepted event, 
both the cross section of Eq. (|13D and the kinematical parameters of the outgoing particles 
are calculated. In Fig. ^ we show the generated distributions of events plotted with respect 
to the invariants of Eq. ([TT|) . 

B. Calculation of energies and scattering angles of the outgoing particles 

Using the definitions of Eq. (^) for the invariants, one obtains the following expressions 
for the energies and scattering angles of the outgoing particles with respect to the beam axis 
in terms of the lab variables: 



2 ^2 



Ek' 



2ml "~ 1 



2me 



E„i 



Si + q — r — mil 
1' ~ 7^ ' 



and 

cos(6k' 



2m,f. 



S2- s - q^ + ml + 2EiEk' 



mu 



'k' 



r'-2ml + 2E,Ef 
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cosi^y^'j 



q^ -r"^ -S2 + ml + 2EiEg 

2\Pi\Eg, 



The distributions generated for the energies and scattering angles are shown in Figs. |] and 
1^, respectively. 

C. Components of the 3-momenta of the outgoing particles 

We define the {x, z) plane as the horizontal plane with the beam axis in the positive z 
direction and the y axis directed vertically, such that Cx, iy, and iz form a positively oriented 
Cartesian basis. The spherical polar components of a generic 3-momentum p are defined as 
usual: 

Px = |p|sin(6')cos(0), 
Py = \p\ sin(^) sin(0), 
Pz = |p|cos(6'). (37) 

The Pz components of the 3-momenta of the outgoing particles are calculated from the ener- 
gies and scattering angles of the outgoing particles. To calculate the px and Py components, 
it is necessary to know the azimuthal angles which are obtained as follows. 



1. From Eq. {^, 

s, = ml + 2Ek'E,, - 2kWx - 2k'X - ^Kq'z, (38) 

we obtain for the azimuthal angles of the final electron and photon 

1 + tan(0fc/) tan(6'g/) cos(0fc/ + 4>q')- (39) 



2Ek'Eg, +ml- si 



2. Analogously, from Eq. ([TT|), 

S2 = ml + 2EfE,. - 2pfX - ^Pfy^l'y " 2^/,?^, (40) 

we obtain for the azimuthal angles of the final pion and photon 
2Ep^E,,+m^-S2 ^ ^ ^ ^^^ ^ 

We generate cf)q' randomly, and then use Eqs. ( pQ]) and (|iT| ) to determine (f)p^ and (f)k'- 
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VII. SENSITIVITY 

In Ref. [Q and in the associated bibliography, detailed descriptions are given of the 
different theoretical models of the pion polarizabilities. Here, we only collect the various 
predictions for the polarizabilities (see Table H). 

For values of a of the order of 10~^^ cm^ (see Table ^ one obtains for the ratio 
Q;/(|7rr^) ~ 10~^, where r^ is the pion electromagnetic radius. A comparison with the 
hydrogen atom, for which this ratio is of the order of 1, leads to the conclusion that the 
pion is a very rigid object, with strong interactions between its charged constituents. As a 
consequence, we expect the effect of the polarizabilities on the RCS or VCS cross sections to 
be rather small. Therefore a crucial component of the event generation is the determination 
of the VCS phase space region, where the data are expected to be sensitive to the pion 
polarizability. 

To that end, we introduce the variable Ratio as 






Here, |A^P is the total squared amplitude, whereas both Re{M.%^j.^M.^-B} and 
Re{A4^^J^Bu} depend linearly and |7WnbP quadratically on the generahzed dipole polar- 
izability aL(g^) [see Eqs. (g), ([MD, (|A5|), and ([A5D]. We generated VCS events for different 
values of a in Eq. (|]). The corresponding distributions of these events with respect to the 
variable Ratio are shown in Fig. |[ For a = 0, we find Ratio values less than 0.05, due to the 
/(g^) term aiiq"^) [see the definition of Eq. (|)]. For a = 2.7 x 10~^^ cm^ and a = 6.8 x 10~^^ 
cm^ the "sensitivity" reaches 20 % and 50 %, respectively. The g^-dependent term /(g^) in 
ttL(g^) decreases relatively with increasing a, as shown in Fig. |^. 
When 

\Ratio\ > 0.05, (43) 

the pion polarizability component of the VCS cross section exceeds 5 %. The simulation 
shows that this occurs only for ~0.1 %, ~0.6 %, and ~1.3 % of the generated events for 
a = 0, a = 2.7, and a = 6.8 x 10^^^ cm^, respectively. 

The differential acceptance of these polarizability-sensitive events is shown as a function 
of invariants in Fig. |T^, from which one finds that the sensitivity of VCS to pion structure 
is significant in the following cases: 

1. the invariant mass Si is larger than 0.1 GeV^, 

2. the absolute value of the squared momentum transfer r^ is larger than 0.05 GeV^, 

3. the invariant mass S2 is in the range 0.05 GeV^ < S2 < 0.2 GeV^. 

Fig. (bottom left) shows that the differential acceptance function versus g^ is approxi- 
mately flat. Therefore, the dependence of the sensitivity to the electron energy is small. 
Fig. (bottom right) shows that the sensitivity increases with increasing absolute value of 
the squared momentum transfer r^. Therefore, the sensitivity is larger for lower values of 
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the outgoing pion momentuin [see Fig. [TT| (top)]. As a result of the infrared divergence of 
the BH amphtude, the sensitivity increases with increasing energy of the outgoing photon 
[see Fig. |1T] (bottom)]. 

Events with Ratio > 0.05 are concentrated in the forward region of the angle 6* in the 
67 CM frame and in the backward region of the angle 6* of the tc'j CM frame, as shown in 



Fig. 12 



As shown by the calculations, by implementing the requirements 



Ef < 0.62Ei, 

Eg, > 0.15Ei, (44) 



one achieves 



SENS = ^ycsia = ^-7)- -ycsia = 6.8) ^ 

avcs(a = 2.7) ^ ^ 

as shown in Fig. [1^ (b). 

However, even when there is sensitivity to a, the sensitivity to the g^ dependence is very 
low. Fig. |1^ shows the difference in the total (sensitive) cross section with and without q^ 
dependence of the A^nb amplitude Eq. (|25| ) as a function of the beam momentum. The 
percentage difference is only of the order of 0.3 %. 



VIII. VCS TOTAL CROSS SECTION 

We have integrated the cross section of Eq. (|T^ for the kinematic region specified in 
Eq. (|3^ ) using two different Monte Carlo integration programs, namely, DIV0N4 [^ and 



VEGAS [^]. The results are 106.7±1.6 nbarn and 105.5±0.26 nbarn, respectively. For 
further calculations, we use VEGAS, since the calculation errors given by VEGAS are about 
six times smaller and the calculation is twice as fast as for DIV0N4. Using VEGAS, we 
calculate the Born cross sections of the reaction Eq. (|I]) and the total cross section of the 
process Eq. (P for various values of a including the g^ dependence given by Eq. (pSj) . The 
results are shown in Fig. ^, and the values are collected in Table fT^. 

The small increase of the VCS total cross section for a = 0, with respect to the VCS- 
Born cross section, is due to the g^ dependent component of the A^nb amplitude [see Eq. 
(p5|)]. The small decrease in the VCS total cross section with increasing a is due to the a 
dependence of the A^nb amplitude [see Eq. (p5|)]. 

We also calculated the vr VCS cross section for a second kinematic region with 

- 0.2 GeV^ <q^ < -0.032 GeV^ (46) 

for the q^ range and ranges for si, S2, and r^ as given by Eq. (^4]). This choice corresponds 
to the data region studied by SELEX for Hadron-Electron elastic scattering [4^. In this 



region, using the VEGAS code, we find 34.7±0.1 nbarn for the vr VCS total cross section. 
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IX. EXPERIMENTAL STATUS 

Pion Virtual Compton Scattering (VCS) via the reaction ire -^ Tve'-f was observed in 
the Fermilab E781 SELEX experiment |^3|. SELEX used a 600 GeV tt~ beam incident on 
target atomic electrons, detecting the incident tt" and the final state tt~ , electron and 7. 
The number of reconstructed events and their distribution with respect to the kinematic 
variables (for the kinematic region studied) were shown to be in reasonable accord with the 
predictions. The limited statistics of the experiment did not allow deducing a. Ref. P3[ 



discussed methods of improving the apparatus for a second-generation pion VCS experiment. 

X. CONCLUSION 

We have given explicit expressions for the differential cross section of the reaction 7r~e — > 
7r~e7, for which the invariant amplitude has been calculated at the one-loop level, 0{p'^), 
in chiral perturbation theory. The non-Born part of the VCS amplitude depends on the 
generalized dipole polarizability ai^q^). In principle, there exist three generalized dipole 
polarizabilities which, however, are degenerate at the one-loop level. The total cross section 
has been calculated using Monte Carlo integration programs. With the VCS event generation 
we were able to identify regions of phase space which are sensitive to the pion polarizability. 
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APPENDIX A: ANALYTICAL EXPRESSIONS OF THE SQUARED INVARIANT 

AMPLITUDES 

In this appendix we list the analytical results for the squared invariant amplitudes as 
obtained after contraction of the relevant tensors. The scalar products appearing naturally 
after contraction can be written in terms of the invariants of Eq. ( pH]) as 

K^ = {k + k'f = Ami - q\ 
P"^ = {Pi+Pff = 4m^-r^ 

/ '- ( 2 2 
? ■ ? = 2 i^ ~^ 

P ■ q = P ■ q' = S2 - ml - -{q^ - r^) 

K-q={k + k')-{k-k') = 0, 

K -q' = Si-ml + -{q^ - r^). 
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K-P = 2s-Si-S2-ml-ml + -{q^ + r^) 

2pi-K = 2s- 2ml - S2 + q^ - ml 
2pf ■ K = 2s - 2si - ml + r^ - S2 

2 / l\2 2 2,2 2 

U2 - m^ = [Pi - q j - m^ = m^ + q - S2 - r . 
The squared contributions read: 



(Al) 



IMbhP 



e^Ffr^l 



n2 r 



{K-q'y-{q-q'y 



xi{K- q')\P ■ q'f + (g ■ g')'[^V + {K ■ Pf] + 2(ir ■ P){P ■ q'){K ■ q'){q ■ q' 
, {K . q'f - (g ■ g')' 



4m2 



\M 



Born I 



e^F{q^ 



n2 



{P^ pg2 + (ir ■ g')' + {q ■ q'f 



1 (ml[q^P^ + S2-U2) + {2p,-Kf 



_l_ J. 



(A2) 



mt 



[S2 - "^^ 



^ ml [q' (P^ - g2 + ^2) + {2pf ■ JT)^] ^ 2^2 ^ ^2 
(2p,-ir)(2p^-ir) + PV^o^2 „2^ 



(m2 - "^^) 



(S2 - ml) {u2 - ml) 



2mi - q' 



(A3) 



2Pe{A^BornA^NB} 



57r?Tl 



^Pe{ai(g2)}P(g2 



m: 



X Upi ■ /sTp/ ■ K + g^m^ 



[g^ — r 



212 



[s2 -ml){u2 -ml 



, U2-ml-q^ r^S2-ml-q^ 

+Pi ■ Kq ■ K — ^ —^ VfKq ■ K- 



S2 - m^ 
-K^q^ - r^) - iqr] 



U2 -m; 



(A4) 



lA^NBp 



e^\ 1 /87rm, 



q^ I Ami ^ ^ 



- |a,(g^)n(g-g')^i^^ + g^(i^-g')' + 2g^(g-g')']. (A5) 



For the evaluation of the Bethe-Heitler-VCS interference term we quote an intermediate 
result: 

4 (^77)2^^^77)2 {^P^.(^ -^K.q' + P.Kq. q') - K,P,{K ■ q')^ 



9 pi' T int 



-Kpq^P ■KK-q' + P.K^q ■ q\q' - q ■ q') - P.q^K ■ q\q' - q ■ q' 
+q'K,P ■ q'q ■ q' - q'P^K ■ q'q' - g' g^P ■ i^g ■ g' 



+9pAiq'-q-q')P-q'K-q' + p-Kiq-q 



/\2i 



(A6) 
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where we have rewritten k = {K + g)/2 and k' = {K — q)/2, and where terms proportional 
to Qp and q' have been dropped since they are contracted with the conserved VCS tensor. 
The interference between the VCS Born terms and the BH terms is given by 



06 



1 



qr^ mj\[K ■ q') -[q-q')\ 



1 



1 



S2-ml U2-ml^ 



A + 



1 



S2-ml U2-ml^ 



B + C 



where 



(A7) 



A = P-Kq- q'[{P ■ Kf + P'^q^] 

-K ■ q'{{P' + q^){K ■q'P-K + q^P ■ q') + P ■ q'[{K ■ q'f - 2q\ ■ q']}, 
B = -K- q'{{q^ - q ■ q')[{P ■ Kf + P'^q^] + P''{K ■ q'f + P-q'P-KK- q'}, 
C = -2[K ■ q'P ■ q\2q^ - q ■ q' + K'^) + P ■ Kq ■ q\q^ + q ■ q' + K'^) ~ P ■ K{K ■ q'f]. 



In Eq. (|A7|) the upper and lower signs refer to a vr"*" and a vr" beam, respectively. 

Finally, the interference between the non-Born terms of the VCS amplitude and the BH 
terms reads 



1 



2Re{M*^^MBu} = ±87rm.i?e{«^(g^)}^F(r^) ^^^^^ ^,^^ _ ^^ ^^^^^ 



x\P-K[{K-q 



'\2/„2 



q-q) + {q-(^Y + {K^ + q'){q-q 



l\2^ 



+P ■ q'K ■ q'[{K ■ q'f + [q^ - q ■ q'W + q ' q') + K\ ■ q'\ 



(A8) 



where the upper and lower signs refer to a tt^ and a n beam, respectively. Here we made 
use of the approximation of Eq. (p5D which is valid at O(p^) in ChPT. 
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TABLES 



Reaction 


a 


Experiment 


7P -^ 77r+n 
77 -^ vr+vr" 


6.8 ±1.4 ±1.2 
20 ±12 
2.2 ±1.6 


Serpukov |p 
Lebedev [§| 
MARK-2 § 



TABLE L Experimental values of a in units of 10 cm' 



Theoretical model 


a 


/9 


Chiral quark model ||35| 


8.0 


-7.8 


Nonrelativistic quark model [36l 


0.05 


- 


NJL model Q 


10.5-12.5 


-(10.3-11.8) 


ChPT [0{p^)] iH 


2.68 ±0.42 


-(2.68 -2.61) ±0.42 


ChPT [0(/)] |22] 


2.4 ±0.5 


-2.1 ±0.5 


Quark confinement model [|38| 


3.63 


-3.41 


QCD sum rules ^ 


5.6 ±0.5 


- 



TABLE II. Theoretical predictions of the polarizabilities in units of 10 ^^ cm^. 





Total cross section (nb) 


Error (nb) 


Born 


105.2 


0.25 


Q = 0.0 


105.9 


0.26 


a=2.7 


105.5 


0.26 


a=6.8 


105.3 


0.25 



TABLE III. Born and total cross sections of the reaction Eq. (|l]) for a negative pion as a 
function of a {a is given in units of 10~^^ cm^). 
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FIG. 1. Three pion Compton scattering reactions: (a) vr Z 



vr Z7, (b) 7p — > 77r'^n, (c) 
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FIG. 2. 0{p^) prediction for the generalized dipole polarizabihties aL{—Q'^) of the charged 
pion (soUd curve) and the neutral pion (dashed curve) as functions of Q^ [see Eqs. @ and (^)]. 
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FIG. 3. Kinematics of TT{pi) + e{k) — > TT{pf) + e{k') + j{q') in the laboratory (lab) frame. 




+ 



FIG. 4. The reaction n {pi) + e^{k) -^ tt (pf) + e^(fc') + "y{q') to lowest order in the electro- 
magnetic coupling: Bethe-Heitler diagrams (a) and (b), VCS diagram (c). The four- momentum 
of the virtual photon emitted by the pion in the BH process is r = pi — pf. The four- momentum 
of the virtual photon absorbed in the pion VCS diagram is q = k — k'. Due to four-momentum 
conservation r = a' — a. 
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FIG. 5. Generated distributions of events plotted with respect to the invariants of Eq. (|llD 
The generation was done with a vr" beam momentum of 650 GeV/c and a = 6.8 x 10~^^ cm^. 
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FIG. 6. Generated distributions of the energies of the final-state particles. The generation was 
done with a 7r~ beam momentum of 650 GeV/c and a = 6.8 x 10~^^ cm'^. 
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FIG. 7. Generated distributions of the 9 scattering angles of the final-state particles. The 
generation was done with a vr" beam momentum of 650 GeV/c and a = 6.8 x 10~^^ cm^. 
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FIG. 8. The distributions of the variable Ratio [see Eq. (l42|)1 for different values of the pion 
polarizability. 
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FIG. 9. Plot of f{q^)/{a) (Eq. (0)) versus g^ for a = 2.7 x IQ-^^ ^nd a = 6.8 x 10 
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FIG. 10. The differential acceptance function of events that satisfy the condition of Eq. (43) 
for the invariants for a = 6.8 x 10~^^ cm^. 



28 



^4 


- 










No cuts 




I 


Ra;io>0.05 


0^ 


- 


1 


0^ 


- 


^^'---^-r- 





r 




p 






1 


1 


r 




'l 




1 1 1 1 1 1 1 1 


( 


} 


1 


c 







200 300 




400 
7T energy 

No cuts 
Ratio>0.05 



)00 



600 



102^ 



10 - 



1 - 



100 



200 



300 



400 



bOO 



600 



7 energy 



FIG. 11. The generated distribution of the energy of the pion (top) and the photon (bottom). 
The dashed line shows the part of the distribution with Ratio > 0.05 [see Eq. (|4^)]. 
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FIG. 12. The generated distribution of the angle 9*^ in the ttj (top) and 67 (bottom) CM 
frames, respectively. Distributions with Ratio > 0.05 [see Eq. ([43|)] are shown on the right. 
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FIG. 13. (a) Total cross section of pion VCS for a = 2.7 x 10~ cm^ as a function of the beam 
momentum including the requirements of Eq. (^4[); (b) SENS distribution as a function of the 
beam momentum including the requirements of Eq. (03) . 
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FIG. 14. Difference in the total (sensitive) cross section as a function of the beam momentum 
with and without q^ dependence of A^nb amphtude [see Eq. (|25|)] for a = 2.7 and a = 6.8 x 10~^^ 



cm^. 
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